Roton-like mode in solid 4 He 
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Recent theoretical and experimental work on solid 4 He has focused on vacancies, or other lat- 
tice defects and whether these can form a Bose-Einstein condensate within the solid. The earliest 
suggestion that this could happen was based on the conjecture that the ground state of the solid 
at T = K would include vacancies. Although no direct evidence for zero point vacancies has 
been presented, a variety of experimental observations have been interpreted as being due to ther- 
mally activated vacancies. Here we directly probe the excitation spectrum of hep solid 4 He using 
inelastic neutron scattering. We identify a branch of delocalized excitations, with both longitudinal 
and transverse dispersion that coexist with acoustic phonons. While the energy gap is larger and 
the characteristic wavevector is shifted to coincide with the position of the (100) Bragg peak, the 
longitudinal effective mass of this mode is within 15% of that for rotons in superfluid 4 He. The 
dispersion relation intersects the longitudinal acoustic phonon half way to the zone boundary, and 
this distinguishes the mode from a conventional optic phonon. Our results provide direct evidence 
for coherent delocalized non-phonon excitations in solid helium, which we identify as delocalized 
vacancy modes. 

PACS numbers: 67.80.B-, 67.75. dt, 61.05.F- 



INTRODUCTION 

Theoretical speculation about possible Bose-Einstein 
condensation and superfluidity in solid He began almost 
40 years ago [H, @, H, 3 . Early experimental searches, us- 
ing macroscopic techniques, found no evidence for the 
'supersolid' state However, evidence for some type 

of phase transition in impure or strained crystals was 
found using ultrasonic propagation 0, Q. Recent ob- 
servations of a decrease in the moment of inertia of the 
solid, in impure or strained crystals, at temperatures be- 
low about 200 mK were interpreted as evidence for a 
supersolid 

El El EH- However, an attempt to force di- 
rect flow through small channels was unsuccessful [13], 
and deep inelastic neutron scattering finds no evidence 
for a zero-momentum condensate fraction 



As in the case of superfluid helium, the excitations in 
the solid could provide microscopic evidence for a pos- 
sible Bose condensate or supersolid behavior if present. 
Indirect measurements of the excitations have been made 
using NMR, ultrasound, pressure, and x-ray and neu- 
tron diffraction. The results are generally inter pre ted in 
terms of vacancy excitation 0, M, [13, 0, El Hj. Va- 
cancies are thermally activated and the activation energy 
increases roughly proportionally to the melting temper- 
ature over a wide range of densities [3, El- However, 
the specific heat inferred from these activation energies 
is greater than the total measured specific heat if they 
are assumed to be localized, and is also too large if the 
vacancies behave like an ideal Bose gas [2l[. A possible 
small contribution of these excitations to the pressure has 
been found [i3| ■ There is no direct experimental evidence 
for of zero-point vacancies in solid helium, and quantum 



Monte-Carlo calculations of perfect crystals do not find 
such vacancies 22 , 23| . 

Here we report inelastic neutron scattering measure- 
ments on low density solid 4 He crystals. The measure- 
ments reveal portions of the spectrum that were not pre- 
viously observable and provide higher resolution than 
earlier experiments [24[ ■ We find no evidence for a phase 
transition within the solid state but we do provide evi- 
dence for dispersive non-phonon excitations. We identify 
these excitations as waves of delocalized vacancies and 
estimate their contribution to the specific heat close to 1 
K. Certain similarities between parts of the vacancy wave 
excitations and the roton spectrum in the superfluid are 
suggestive of strong links between the solid and liquid 
phases. 



METHODS 

The excitation spectra were measured using the NG4 
disk chopper time-of-flight neutron spectrometer DCS at 
the NIST Center for Neutron Research, Gaithersburg, 
MD j2f|. The DCS uses choppers to create pulses of 
monochromatic neutrons. Energy transfers on scatter- 
ing are determined from arrival times in detectors that 
range in scattering angle from -30° to 140° in the hori- 
zontal plane and subtend a solid angle of « 0.71 sr at the 
sample. Unless otherwise mentioned, we operated in the 
low resolution setup [ID] which gives a resolution of w 0.2 
meV full-width half-maximum (FWHM) at an incident 
wavelength of using 3.7 A wavelength (energy 6.0 meV). 

The 4 He crystals were grown in a stainless steel cylin- 
drical sample cell with the axis along the vertical. The 
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inside diameter of the cell was 20.5 mm and the outside 
diameter was 22 mm. The neutron beam was restricted 
to a height of 25 mm (smaller than the sample height) 
using a cadmium mask. The total volume of the cell was 
approximately 30 cm 3 . The lower flange of the sample 
chamber was machined out of Cu and included a ring of 
sintered Cu powder to enhance thermal contact to the He 
sample. Contact between this Cu flange and the mixing 
chamber of the dilution refrigerator was made with an an- 
nealed copper ribbon 1.9 mm wide by 0.55 mm thick. The 
upper flange was made of brass and provided a mechan- 
ical mount to the mixing chamber. X-cut quartz trans- 
ducers (resonant frequency 9.77 MHz) were mounted at 
the centers of the top and bottom flanges to measure 
the acoustic velocity and attenuation through the helium 
along the vertical axis. Two additional transducers were 
mounted in the top flange outside the neutron beam, to 
measure the acoustic properties in a horizontal direction. 
The helium prepared for the samples contained w 0.004 
% 3 He (40 ppm) at the start of the experiment but was 
diluted by a factor of 5 to 10 by the time the neutron 
measurements were made. 

Temperature was measured by a R11O2 resistance ther- 
mometer attached to the mixing chamber. The base tem- 
perature reached by that thermometer was 55 mK. The 
ultrasonics were used to monitor crystal growth, as an 
internal thermometer of the solid helium, and as an in- 
dicator of crystal quality. The large size of this sample 
relative (w 30 cm 3 ) to that of previous acoustic mea- 
surements (« 1 cm 3 e.g. Ref. 0) and the long acoustic 
path (44.5 mm versus 8.43 mm) made it difficult to grow 
perfect crystals through the entire acoustic path. In per- 
fect crystals, the temperature dependence of the velocity 
is given by " (T J~^° K) = - 6.3 xl0~ 4 T 4 [ljj, and the 

acoustic attenuation, a, varies as a a T 3 / 2 exp(— A/T), 
where A is associated with the vacancy activation en- 
ergy. In strained crystals, the attenuation shows a broad 
maximum around 0.5 K and the fractional change in the 
velocity from T<ClKtoT = lKis much more than 
6.3 xl0~ 4 and departs strongly from the T 4 dependence. 
In the first crystal measured in this work, the fractional 
velocity change between 1 K and base temperature was 
3.4xl0~ 3 , indicating significant strain. The T depen- 
dences of v and a were also characteristic of crystals 
containing non-zero dislocation densities Q ■ This crystal 
also displayed a small acoustic anomaly like that found in 
Ref.0, confirming that measurements were made at tem- 
peratures below the acoustic anomaly at 200 mK in solid 
hep 4 He. The crystals examined in the work described 
here were grown from the liquid at nearly constant pres- 
sure and temperature, not by the blocked capillary tech- 
nique. 

The first sample (Crystal 1) measured was grown with 
the mixing chamber at 1 K. The sample grew as a sin- 
gle crystal with the [001] axis vertical to within < 1°. 



The a lattice parameter was 3.68 A. The absolute value 
of the acoustic velocity (450 ms _1 ) was also consistent 
with this orientation [26|. Using the spectrometer in a 
white-beam (Laue) mode all peaks observed were associ- 
ated with the c-axis oriented single crystal. The second 
sample was grown with the mixing chamber at 0.5 K 
and consisted of several large crystals, oriented in differ- 
ent directions. The orientation of one of these crystals 
(Crystal 2) permitted acquisition of clean spectra of the 
longitudinal excitations about the (100) Bragg peak. The 
temperature inside the sample chamber during growth of 
the crystals was higher than that of the mixing chamber 
since the heat leak through superfluid in the filling capil- 
lary is much higher than when the capillary is filled with 
solid helium. Measurements with a thermometer inside 
a similar sample chamber indicate that the growth tem- 
peratures were approximately 1.2 K. In both cases, the 
observed spectra showed no evidence of liquid helium in 
the path of the beam. Spectra were also obtained for the 
superfluid at 24 bar, the empty cell at T = 5 K and a 
vanadium standard was used to calibrate detector effi- 
ciencies. Unless otherwise stated all errors refer to one 
standard deviation. 



RESULTS 

Figures QJt and |T|b show scattering intensity contours 
in energy-wavevector space for wavevectors in the basal 
plane, both transverse and longitudinal relative to the 
(100) Bragg peak for Crystal 1. The a lattice parameter 
was 3.69 A. A range of incident energies were used, to 
provide greater coverage in momentum space or better 
energy resolution. For each incident energy, at least 21 
different orientations of the crystal relative to the inci- 
dent beam were used, covering a rotation of 42° about 
the axis perpendicular to the scattering plane. For com- 
parison Fig. [Tfc shows the corresponding inelastic spectra 
for the superfluid at a pressure just below the freezing 
pressure (24 bar). The dominant new features in the 
longitudinal cut (Fig. QJi) are: 

• a quadratic dispersion curve centered at the (100) 
reciprocal lattice point. 

• the apparent intersection and subsequent disap- 
pearance of the two sharp longitudinal modes about 
half way to the zone boundary. 

• a marked asymmetry in the intensity of the acoustic 
phonon which reduces the point group symmetry to 
triangular (24. 27|. 



• a broad spectral feature that appears at about 2 
meV, for wavevectors beyond the point of inter- 
section mentioned above, and extends beyond the 
energy and wavevector range measured here. 



3 



> 
n 

E 



3 
2 
1 

°0.0 



0.5 1 .0 

[h,0,0] 



1.5 



E 



-0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 
[1-h,2h,0] 



^C) 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 
IQI fA' 1 ) 



40 



30 



1.5 



u 



20 



10 



40 



30 



20 



I 



10 



40 



30 



20 



10 



Figure 1: (a) Cut along the [h,0,0] direction, integrating over 
the transverse direction by ± 0.01 reciprocal lattice units 
(r.l.u.). (b) Cut of points (l-h,2hfi) (integrating over ± 0.01 
r.l.u.). These data were taken at base temperature, with an 
incident wavelength of 3.7 A. (c) Superfluid data, taken at a 
pressure of 24 bar and ~ 0.22 K with an incident wavelength 
of 3.49 A. All data have been treated using empty cell and 
vanadium data as described. The (100) Bragg peak has an 
absolute momentum transfer of 1.97 A -1 . 
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Figure 2: Neutron weighted density of states integrating over 
all measured momentum transfers at base temperature (closed 
circles) and high temperature (open circles) for an incident 
neutron wavelength of 5.2 A. Using the thermometer placed 
on the mixing chamber, base temperature is 0.055 K and high 
temperature is 0.6 K. 



momentum transfers. Comparing Figs. [TJt and [TJ:, the 
spectra appear superficially to be very similar. Clyde 
and Griffin [29[ have argued that, in the superfluid, the 
maxon and roton can be interpreted as single particle ex- 
citations, closely coupled to the existence of a superfluid 
condensate. 

No differences in the excitation spectra are observed 
(Fig. [2]) above and below the temperature where both an 
acoustic anomaly and the onset of a decrease in moment 
of inertia were observed (« 200 mK). The longitudinal 
data were all confirmed on Crystal 2. 

Some high energy resolution data was also taken to 
look for (dispersionless) excitations with energies « 1 K. 
Such modes had been found using ultrasound We 
could not resolve any such excitations in measurements 
spanning approximately half the first Brillouin zone and 
the energy range from meV to 1.2 meV. The nominal 
energy resolution for the spectrometer in this configura- 
tion was 54 /icV (FWHM) for an energy transfer of 
meV. We note that in Ref. LUJ this low activation energy 
was only seen in crystals grown below 1 K. The growth 
temperature in the current work is estimated to be 1.2 K 
or slightly higher. 



This latter broad feature is also seen in the superfluid 
(Fig. [TJ:) , starting at a similar absolute momentum trans- 
fer (ss 2.4 A -1 ) and a slightly lower energy, corresponding 
to approximately twice the minimum energy gap in each 
case. This feature is also seen in the normal fluid 28[. In 
the solid, the high energy excitation displays hexagonal 
symmetry about the origin. In Ref. [28|, this scattering 
was identified as broad, multiphonon scattering. Alter- 
natively, it might also be associated with strongly inter- 
acting single-particle scattering, in which case it should 
connect with the free particle recoil scattering at larger 



VACANCY EXCITATIONS 

We now examine the low energy excitations close to 
the (100) Bragg point more carefully Figure [3] shows 
the dispersions of these excitations, as obtained by fit- 
ting Gaussians to the data. The lowest energy mode 
is an acoustic phonon and we argue that the sharp in- 
termediate energy mode is a collective vacancy excita- 
tion. Minkiewicz et al. 24( associated the intermediate 



energy mode with an optical phonon despite the dis- 
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Figure 3: Excitations in the longitudinal and transverse directions from the (100) Bragg peak, measured at 0.055 K with a 
neutron wavelength of 3.7 A for the longitudinal cuts (a,b) and 5.2 A for the transverse cuts (c,d). The positions (a,c) and 
intensities (b,d) were obtained from Gaussian fitting. The acoustic phonon branches are shown by filled circles, the vacancy 
mode by open circles, the optic branches by open squares, and in the transverse case there is an additional mode marked by 
open triangles. The neutron wavelength was 3.7 A for the longitudinal cuts (a,b) and 5.2 A for the transverse cuts (c,d), apart 
from the optic mode in the transverse direction, which was measured at 3.7 A. For comparison the vacancy mode at 3.7 A is also 
shown (asterisks) in (b) and used to scale the intensity of the 3.7 Adata in (d) to that at 5.2 A. In the upper panels (a,c), the 
lines are the calculated dispersion curves for the acoustic (dashed) and optic (dotted) phonon branches from the self-consistent 
phonon calculations in Ref. [3(| In (b) , the intensities are calculated from the inelastic structure factor following Ref. [3(] with 
an additional common scaling factor. 



crcpancy with the self-consistent phonon calculations of 
Gillis, Koehler and Werthamer |30[, to be discussed be- 
low. The intensities for the longitudinal branches were 
calculated from the inelastic structure factor based on 
this self-consistent phonon theory pTl |30T | , and compare 
favorably to the data in Fig.[3]except at the intersection. 
In particular, the asymmetry of the acoustic phonon 
branch is well reproduced. 

We present two arguments for associating the interme- 
diate energy mode with a vacancy excitation, (i) This 
mode intersects with the longitudinal acoustic phonon 
branch. On symmetry grounds, longitudinal optic and 
acoustic phonon branches in a hep structure may not 
cross or be degenerate except at the zone boundary in 
the c* direction 31| . For the transverse modes, where 
the zone boundary occurs at h — 0.33, there is no col- 
lision between the strongest modes, presumably due to 
the lower transverse sound velocity. 

(ii) If the positions of the optic and acoustic branches 
are taken from the self-consistent phonon calculations 
and the existence of a third dispersive excitation is as- 
sumed, excellent fits to the data can be obtained using 
damped harmonic oscillator lincshapes convoluted with 
a Gaussian instrumental resolution function of 0.2 meV 
FWHM (Fig. 0J. The optic branch corresponds to the 



broad, weak excitation at higher energies, with a signifi- 
cantly broadened lincwidth. This assignment resolves the 
discrepancy between the self-consistent phonon calcula- 
tions and experimental observations. Close to their inter- 
section both the acoustic branch and the vacancy mode 
widen significantly, indicating decay of phonons into va- 
cancy waves. Above the intersection, just one sharp mode 
persists. 

Markovich et al. [32| observed an optic-like mode in the 
bec phase of 4 Hc in a longitudinal cut through the (110) 
Bragg peak, although none should exist in a monatomic 
bec structure. The strong similarity between this spec- 
trum and that reported here in the hep phase suggests 
that this branch of the spectrum in the bec phase could 
be the vacancy excitation. 

With our identification of the intermediate energy 
mode as a vacancy excitation, the minimum energy is 
identified as the vacancy activation energy, and is found 
to be 11 K. This is consistent with values obtained by 
NMR [lj| and observation of activated thermal expan- 
sion above ~ 10 K in x-ray 17[ and neutron diffraction 
flit studies. The higher activation energy observed in 
the bec phase [l6| is consistent with the higher melting 
temperature of the bec crystal. The dispersion of this va- 
cancy mode is quadratic, and so an effective mass can be 



5 




12 3 

Energy transfer (meV) 



Figure 4: Energy scans at (hOO) integrating over a momen- 
tum transfer of h ± 0.005 r.l.u. in the longitudinal direction 
and ± 0.01 r.l.u. in the transverse direction, for three selected 
values of H. The wavelength used was 3.7 A, at base tem- 
perature. The solid lines are fits using a model consisting of 
three damped harmonic oscillators convoluted with a Gaus- 
sian resolution function of 0.2 meV FWHM. The positions of 
two of the oscillators are fixed using the phonon dispersions 
calculated by Gillis, Koehler and Werthamer [3C| . and the 
third (middle) oscillator determines the position of the va- 
cancy mode. At h = 1.18 the lower and intermediate energy 
peaks have merged into a single peak and the spectrum can 
be described by two modes. 
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Figure 5: Energy cuts at (l-h,2h,0) integrating over a momen- 
tum transfer of h ± 0.005 r.l.u. in the transverse direction and 
± 0.01 r.l.u. in the longitudinal direction, for three selected 
values of h. The lines are Gaussian fits to the data. The 
wavelength was 5.2 A, and these measurements were taken 
at 650 mK. At h = -0.11 the highest energy mode is out of 
the accessible energy range. As h decreases, the intensity of 
the middle mode decreases, finally disappearing at h = -0.07. 
The spectra are displaced vertically by 35 units. 



extracted. Parameters for this quadratic dispersion are 
given in Table |TJ The uncertainties quoted include only 
those from the fit and therefore may underestimate the 
total uncertainty. This excitation was resolution limited 
and from our highest resolution data with incident neu- 
tron wavelength 5.2 A we can place an upper bound of 
21 /ieV on the linewidth. The possibility of a dispersive 
mode associated with defects (specifically 3 He impurities) 
in 4 He has been discussed by Manousakis [3^|, although 
we note that the amount of 3 He in our sample is less than 
0.0008%. 

To summarize, we have observed a gapped zone-center 
excitation that is not a conventional optic phonon, be- 
cause the longitudinal optic and acoustic modes cannot 
cross and annihilate on symmetry grounds, and because 
a similar mode was observed in the bec phase (32| , which 
has no optic mode. 

An additional mode is observed between the transverse 
acoustic and vacancy modes (Fig. [5]). This new mode 
possesses the lattice symmetry and therefore cannot be- 



long to a second crystal, a possibility also dismissed by 
Laue scans of the crystal. It is also observed at two differ- 
ent incident wavelengths, 5.2 and 3.7 A, indicating that 
it is not due to multiple scattering. It is natural to asso- 
ciate both of these non-phonon features of the hep 4 He 
excitation spectrum with some form of propagating de- 
fects or vacancies. 

In hep 4 He, the specific heat capacity is not fully ex- 
plained by phonon contributions [341 ] . and it seems likely 
that this excess specific heat belongs, at least in part, to 
the delocalized vacancy excitations seen here. Because 
of the substantial gap and small effective mass, thermal 
population of these modes does not account for an appre- 
ciable part of the excess specific heat. It would appear 
that the T 7 term in the measured specific heat is due 
almost entirely to lattice anharmonicity. This problem 
may ultimately be resolved by determination of the com- 
plete density of states for all excitations, whereas so far 
we have only measured the spectra for wavevector trans- 
fer in the basal plane. 
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Mode 


Effective mass 
niHo 


Energy gap 
meV 


Vacancy minimum 

A- 1 


Longitudinal 
Crystal 1 


0.107 ± 0.003 


0.951 ± 0.004 


1.964 ± 0.006 


Longitudinal 
Crystal 2 


0.114 ± 0.005 


1.002 ± 0.003 


1.970 ± 0.006 


Transverse 
Crystal 1 


0.144 ± 0.002 


0.916 ± 0.002 


1.970 ± 0.001 


Roton 
Superfluid (24 bar) 


0.094 ± 0.002 


0.680 ± 0.001 


2.054 ± 0.001 



Table I: Characteristics of the vacancy mode in the solid and the roton mode in the superfluid. 



CONCLUSIONS 

A new coherent excitation has been observed in single 
crystalline 4 He. We argue that it is the elementary exci- 
tation associated with delocalized vacancies in 4 He and 
set a lower limit on its lifetime of 0.2 ns. This new exci- 
tation, as interpreted here, makes an insignificant contri- 
bution to the excess specific heat at high temperatures. 
We have also observed a broad feature in the spectrum 
which appears to be the optical phonon branch that was 
not observed in previous measurements. The longitu- 
dinal acoustic phonon and vacancy modes intersect and 
merge into a single broadened mode, with a linear dis- 
persion at high energies. 

While to high accuracy there is no change in the den- 
sity of states between 0.15 meV and 5.5 meV upon cool- 
ing through the putative supersolid transition, the new 
collective modes that we have observed remain an in- 
triguing aspect of solid helium. The similarities between 
the superfluid spectrum and the longitudinal part of the 
solid spectrum are remarkable and indicate that roton- 
like entities persist in the solid state, modified principally 
by the appearance of transverse dispersion. 
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